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The electromagnetic properties of superconductors near a ferromagnetic instability are investi- 
gated by means of a generalized Ginzburg-Landau theory. It is found that the magnetic flux ex- 
pulsion capability of the superconductor gets stronger, in a well-defined sense, as the normal-state 
magnetic susceptibility increases. The temperature dependencies of the London penetration depth, 
the critical fields, and the critical current are all strongly affected by ferromagnetic fluctuations. In 
particular, for the critical current we find a temperature exponent a ~ 2 over an appreciable tem- 
perature range. The extent to which proximity to magnetic criticality may be a viable explanation 
for recent observations in MgCNi microfibers, which find a ~ 2, is discussed. 
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I. INTRODUCTION 

The coexistence of ferromagnetism and supercon- 
ductivity has received substantial attention over the 
years. Around 1980, such states were predictedi*^ and 
observed^^ and the topic later received renewed inter- 
est in the context of experimental observations in rare 
earth borocarbides£ More recently, interest in this sub- 
ject has been revived by the observation of coexisting 
superconductivity and ferromagnetism in UGe^i^ and 
URhGc, 9 where both types of order are believed to 
be due to electrons in the same band. Recent theo- 
retical attention has centered on the structure of the 
phase diagramjiS on the existence of spontaneous flux 
latticeSfiiiiii£ and on the question of spin-triplet versus 
spin-singlet superconductivity^ 

In contrast, less is known about the properties of su- 
perconductors on the paramagnetic side of, but close to, a 
ferromagnetic instability. We will refer to "paramagnetic 
superconductors" to describe systems in this regime, al- 
though the superconductivity of course leads to the usual 
strong diamagnetic effects. Such paramagnetic super- 
conductors include systems below the superconducting 
transition temperature, but above the temperature be- 
low which coexistence of superconductivity and ferro- 
magnetism occurs, as well as systems that never develop 
ferromagnetism, but are close to a ferromagnetic insta- 
bility in some direction in parameter space other than 
temperature. An example of the latter is believed to be 
the non-oxide perovskite MgCNi3, which superconducts 
below a critical temperature T c w 8 Km^ There is no ev- 
idence for a ferromagnetic phase in this material, but 
it has been suggested that a ferromagnetic ground state 
can be reached upon a relatively small amount of hole 
doping. 15 This system may thus be close to a ferromag- 
netic instability everywhere in its superconducting phase. 

A recent study of MgCNi3 microfibers, with T c = 
7.8 K, has revealed an anomalous temperature depen- 



dence of the critical current density j The critical 
current density vanishes at T c according to a power law 
j c cx \T - T c \ a , with a = 2 between about 1% and 10% 
away from the critical point, and no crossover to the usual 
Ginzburg-Landau behavior, which predicts a — 3/2. The 
authors of Ref . have ruled out morphological effects as 
an explanation, which raises the question whether prox- 
imity to a ferromagnetic state may be responsible. In- 
deed, since ferromagnetic fluctuations are expected to 
weaken (singlet) superconductivity, this is a plausible 
suggestion for the origin of the weaker-than-expected 
temperature dependence of j c . 

The probable proximity to ferromagnetism has led to a 
debate about the nature and symmetry of the pairing in 
MgCNi 3 ^i This point has not been settled; some exper- 
imental evidence points to conventional s-wave pairing; 
other, to a superconducting order parameter with nodes. 
The nature of the pairing in the other materials men- 
tioned above has not been unambiguously determined ei- 
ther. In this paper we will focus on the behavior close to 
T c , which is qualitatively independent of the symmetry 
of the order parameter and thus expected to be the same 
for all nearly ferromagnetic superconductors. We use a 
generalized Ginzburg-Landau theory to theoretically in- 
vestigate the electrodynamic properties of a supercon- 
ductor as a ferromagnetic instability is approached. We 
treat the superconductivity in the usual mean-field ap- 
proximation, but the magnetic critical behavior exactly 
in a scaling sense. Somewhat counter-intuitively, strong 
magnetic fluctuations make, in a well-defined sense, the 
superconductivity more robust in certain respects. In 
particular, the penetration depth becomes anomalously 
short. The thermodynamic critical field, on the other 
hand, becomes weaker, as one might intuitively expect. 
The temperature dependencies of the critical field H c and 
the penetration depth A depend on the magnetic critical 
exponents 6 and 7, respectively. For the critical current 
j c oc H c /X, this results in an exponent a between 1.5 (the 
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Ginzburg-Landau result) and 2.16 in various temperature 
regimes. We will discuss both the existing experimental 
observations, and predictions for new experiments, in the 
light of these results. 

This paper is organized as follows. In Sec. [n] we give 
elementary phenomenological arguments for the depen- 
dence of the thermodynamic critical field, the penetra- 
tion depth, and the critical current density, on a con- 
stant normal-state magnetic permeability fi n . We then 
generalize these results to the magnetically critical case, 
where one needs to distinguish between /j, n and the spin 
susceptibility fi s in a superconduting state, and both [x n 
and fi s become nonanalytic functions of various control 
parameters. In Sec. IIIII we derive these results from a 
generalized Ginzburg-Landau theory, and in Sec. II VI we 
give a discussion of our results. 



II. PHENOMENOLOGICAL ARGUMENTS 
A. Paramagnetic systems 

We are interested in the electromagnetic properties of 
superconductors with ferromagnetic fluctuations. We de- 
note the normal-state spin susceptibility, which describes 
the response of the spin degrees of freedom to an ex- 
ternal magnetic field in the absence of superconductiv- 
ity, by Xni smd the corresponding spin permeability by 
fj, n = 1 + 47rx n . This is in contrast to the spin perme- 
ability fi s = 1 + 47rx s , which includes the effects of the 
superconductivity on the spin response, and the magnetic 
permeability /i = l+Anx, which describes the response of 
the total magnetization, including the diamagnetic part. 
It is instructive to first recall the dependence of super- 
conducting properties on a constant /i„ ^ 1, neglecting 
the distinction between /x n and fj, s m^i^ This can be done 
by means of elementary arguments. 



1. Thermodynamic critical field 

Consider the free energy density / of a system in a 
magnetic field. It obeys 

df = df(H = 0) + -^HdB, (2.1) 

where H is the thermodynamic magnetic field, and B 
is the magnetic induction. For the sake of simplicity, 
we ignore the vector nature of various quantities in our 
free energy considerations. For fixed B, f is the appro- 
priate thermodynamic potential whose minimum deter- 
mines the equilibrium state. However, in an experiment 
H is fixed, since (c/4-7r)V x H = j ext is the external cur- 
rent density, and only the latter is experimentally con- 
trolled. One therefore must perform a Legendre trans- 
form to a thermodynamic potential g = f — BH/Att^SL 



which obeys 

dg = df(H = 0) - — B dH. (2.2) 
Air 

In a paramagnetic phase, including paramagnetic super- 
conductors, the relation between B and H is 

B = H + AnM = (1 + An X )H = fiH, (2.3) 

with M the magnetization, x(^j H) the magnetic sus- 
ceptibility, and fi = 1 + Att\ the magnetic permeability. 
Integration of Eq. 12.211 yields 

1 f H 

g(T, H) = f(T, H = Q)-—J dh [1 + 4ttx(T, h)] h. 

(2.4) 

This is generally valid. In a superconducting Meissner 
state, 5 = 0, and hence \ = — l/Air (ideal diamag- 
netism), and f(T, H = Q) = f + t\ip\ 2 /2 + u\tp\ 4 /A, with 
/o the free energy density of the normal state, ip the 
superconducting order parameter, t cx (T — T c )/T c the 
dimensionless distance from the superconducting critical 
point, and u a parameter. In a normal metal far from a 
magnetic instability, and ignoring normal-state diamag- 
netic effects, x(^) H) ~ const. = Xm or /i n = 1 + 4-7rx n = 
const., and f(T, H = 0) = /o- In a normal metal close to 
a ferromagnetic critical point, x is a complicated function 
of T and H. 

Now consider a superconductor with /Lt n = const. Ac- 
cording to Eq. H2.4fl . the magnetic energy density gained 
by the system allowing magnetic flux to penetrate, i.e., 
the free energy density difference between the Meissner 
state with B — and the normal state with B = /i n H, 
is E m /V — n u H 2 /8tt. By contrast, the condensation en- 
ergy density gained by the system becoming a supercon- 
ductor is E con d/V — t 2 /Au. The thermodynamic critical 
field, which is defined by these two energies being equal, 
is thus 

H c = y^J^i |t|A//^ - H°J^, (2.5) 

with H® = \j2n ju \ t\ the critical field for a system with 
/i n = 1. An increase in /j, n thus decreases the critical field, 
as one might expect since the externally applied field is 
amplified inside the material. 

2. London penetration depth 

The dependence of the London penetration depth A 
on /i n is intuitively less obvious. Consider a large super- 
conducting sample, with linear dimension L, surrounded 
by vacuum and subject to a homogeneous external mag- 
netic field H = (0, 0, H) in z-direction. Along the left 
edge of the sample, the magnetic induction will be of 
the form B(x) = (0,0, B{x)) with B(x) = B Q e~ x ' x 
(x > 0). To determine Bq, imagine a thin (thickness d) 
layer of normal conducting material around the super- 
conductor. Except for the superconductivity, the normal 
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layer should have the same properties as the supercon- 
ductor, in particular, a magnetic permeability fi n . Then 
we have B = H in vacuum, and B = [i n H inside the 
normal layer, see Fig. QJa). Now let d — ► 0. Then we 





d 



FIG. 1: Magnetic induction schematically as a function of 
position at a vacuum (V) - normal metal (N) - superconductor 
interface (S) (a), and at a vacuum - superconductor interface 
(b). 



have (Fig.njb)) 



B(x) 



H for x < 
[i n He- x / x for x>0 



(2.6) 



Now consider the current density associated with B(x). 
From Ampere's law we have 



j(x) = — V x B(x) 
with c the speed of light and 



(0,j(x),0), 



i(*) = i^*(*>- 



(2.7a) 



(2.7b) 



This is the total current density. It has three contribu- 
tions, namely, the supercurrent density j sc , the spin or 
magnetization current j S pm = cV x M , with M the spin 
contribution to the total magnetization, and the external 
current density j ext = cVxiJ/47r. The latter vanishes in 
the case we are considering. In a normal metal, the spin 
current is the only contribution if we ignore normal-state 
diamagnetic effects. The spin or normal-state suscepti- 
bility Xn is defined as the response of M to the total 
magnetic induction B minus the contribution to B of 
M itself, 

M = X n(B - AirM) = ( Xa /» n ) B. (2.8) 
For the supercurrent density j sc = j — j sp i n this implies 



47T^ n 



V x B(x) = (0,j BC (x),0), 



with 



B(x). 



(2.9a) 



(2.9b) 



47T/U n A 

Now consider one surface (area L 2 ) of the sample. Ne- 
glecting corner effects, and for A <C L, the total magnetic 
flux $ through that surface is 

$=/ dy / dxB(x)K,LI dxB(x) = LXfx n H. 
Jo Jo Jo 

(2.10) 



On the other hand, the total supercurrent flowing near 
that surface is, from Eq. I|2.9b|) . 



L, I dxj sc mL j dxj S c(x) 



L 



47T/i n A 



We thus can write the flux 



47TjJ n A 

c L 



47T/in A 

Nqv, 

c L 



(2.11) 



(2.12) 



where N is the number of supercurrent carrying particles, 
q is their charge, and v is their velocity. If m is their mass, 
then Ekm = Nmv 2 /2 is the total kinetic energy of the 
supercurrent. The flux can thus be written 



47T/Z n A 

c L 



-''kin ■ 



(2.13) 



Now we make two observations. First, N ps L 2 Xn, with 
n = \ip\ 2 the particle number density. Second, at the 
critical field strength the kinetic energy of the supercur- 
rent must equal the condensation energy in the region 
where the current is flowing, which is (see Sec. Ill A|) 
-Etond = L 2 Xt 2 /u. With A the London penetration 
depth for /i n = 1, 



A = ^/mc 2 j^q 2 \^\ 2 cx l/\t\ 1/2 , 
this allows to write the flux at the critical field 
^ c = XLH^ a X/Xo = LXfi a H c , 



(2.14) 



(2.15) 



where the first equality follows from Eq. I|2.13|l . and the 
second one from Eq. I|2.10|l . We thus obtain 



,,18.19 



H c /X — -ff,?/Ao , 



A = A /a//m 



(2.16a) 



(2.16b) 



The penetration depth thus decreases with increasing /i n , 
as does the critical field. This is somewhat counterintu- 
itive, as it implies that the superconductivity becomes 
in some sense more robust. It also implies that a large 
normal-state magnetic permeability will make the super- 
conductor necessarily of type Ii£ We will come back to 
this observation. 

Notice that the above derivation relies only on very 
general energetic considerations and on Ampere's law. 
Also notice that it uses an identity at the critical field 
strength, where the superconductivity vanishes. This is 
fine for H c , but the penetration depth is a property of 
the superconducting state, and hence the use of /x n is not 
quite appropriate for this quantity, except in the limit 
A — > oo. More generally, A depends on ji s , which in turn 
depends on the superconducting properties. This makes 
no difference deep inside the paramagnetic superconduct- 
ing phase, and Eq. I|2.16bfl is valid there. However, as we 
will see it makes a crucial difference close to a ferromag- 
netic instability. 
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3. Critical current 

In order to discuss the critical current, we assume a 
thin-wire geometry with wire radius i?^2*Si The super- 
current density, which is the total current density minus 
the spin current density, can be written as a generaliza- 
tion of Eq. (|2.9a[l . 



(2.17) 



j sc (x) = — [V x B(x) - 4ttV x M{x)\ 
V x B(x), 



&TTf/,(x) 



where we have used Eq. I|2.3|) and n(x) is the local mag- 
netic susceptibility. Now integrate over the cross section 
of the wire. Assuming a homogeneous current density 
within a distance A from the surface, and using Gauss's 
theorem on the right-hand side, we have 

27ri?A:/ sc = — — &dt-(V x B(x)) = -RH 

where we have used Eq. I|2.6|) . The critical current den- 
sity j c is the one that produces the thermodynamic crit- 
ical field H c , which yields the familiar London theory 
result 



f 20 



cH c /4tt\. 



(2.18) 



This result is plausible: Dimensionally, j c must be a mag- 
netic field divided by a length. The relevant length scale 
is the thickness of the area that supports diamagnetic 
currents, which is A. The relevant field scale should be 
the field that corresponds to the condensation energy, 
which is H c . To the extent that /j, s w ^ n = const., as we 
have assumed in Sec. Ill A~2l Eq. (|2.1tia[l implies that j c 
is independent of /z n , 



Jc 



Jc 



(2.19) 



As we will see below, this result changes drastically in 
the vicinity of a ferromagnetic instability. 



theory, <5 = 3M Substituting Eq. lf2~2"0|) into Eq. (|2~5|) . 

we find for the thermodynamic critical field 



H c oc |t| 2 <V(*+D 



(2.21) 



This result holds for a system where the distance t from 
the superconducting critical point can be changed while 
the system remains tuned to magnetic criticality (more 
precisely, to the parameter values where magnetic criti- 
cality would occur in the absence of superconductivity). 
Generically, the dimensionless distance r from magnetic 
criticality will change as well if t is changed, and we will 
discuss such more realistic situations in Sec. IIVI 

For the penetration depth, the situation is more com- 
plicated. In contrast to H c , which compares the normal- 
state magnetic energy with the superconducting conden- 
sation energy that has nothing to do with spin mag- 
netism, A is entirely a property of the superconduct- 
ing state, and the feedback of the superconductivity on 
the spin susceptibility, or the difference between /j, n and 
fi a , cannot be neglected. As a result of this feedback, 
the magnetic transition in the presence of superconduc- 
tivity does not occur at r = 0, but rather at a value 
r oc — £m/Aoi Here is the magnetic correlation length 
at zero temperature. This suggests that the spin sus- 
ceptibility at r — will be effectively Xs oc Ao/£m S> 1 
in a mean-field approximation. More generally, one has 
fi a w Xs oc (Ao/i;^) 7 , with 7 another critical exponent. 
Using this in Eq. I|2.16b(l with /u n replaced by fi s , we ob- 
tain 



A oc Aq 



-7/2 ( C \l/ 2 



ie^r 2 oc ir l/2+ ^ /4 . 



(2.22) 



Since 7 fs 1.4 > for ferromagnetic systems^ this im- 
plies that the penetration depth at magnetic criticality is 
anomalously short. Close to the superconducting transi- 
tion, the superconductor will therefore also be of type I, 
in agreement with a conclusion drawn from studying the 
ferromagnetic phased 

For the critical current density, Eqs. (|2.18() . (|2.21() . and 
(|2~2lt predict 



Jc oc 



\t\° 



(2.23a) 



B. Systems at a ferromagnetic instability 



with 



In the vicinity of a ferromagnetic instability of the nor- 
mal metal, the normal state magnetic susceptibility Xm 
and hence the permeability /i n , become large and diverge 
as the phase transition is approached. At a ferromag- 
netic critical point, the region of linear response shrinks 
to zero, and Xn and ^ n become strongly field dependent. 
This field dependence is characterized by the critical ex- 
ponent 5j2£ 

Mn^XnOci/ 1 /*- 1 . (2.20) 

The value of 5 depends on the universality class the par- 
ticular magnetic system belongs to. For all realistic uni- 
versality classes, 5 ~ 5, whereas in Landau or mean-field 



a = 2tf/(£ + l) + l/2-7/4 , (r = 0) . (2.23b) 

With 8 ~ 5 and 7 w 1.4 this yields a 1.8, in contrast 
to the Ginzburg-Landau result a = 3/2. 

These results hold at r = 0, and again we have assumed 
that t can be varied independently of r. Let us relax 
the former condition. From the above argument for the 
effective value of Xs at r = it also follows that Eq. 
(|2.22l) is valid only for |rj < Since ^ is typically 

on the order of a few A, while Ao is typically several 
hundred A or even larger even at zero temperature, and 
diverges as |i| -1 ^ 2 for t — > 0, this is a very small range. By 
contrast, Eq. (|2.2U|) can be valid for r as large as several 
percent, provided H is not too small. Not too close to 
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T c , where H c goes to zero, Eq. (|2.21|) can thus be valid 
in a substantial r-range, while A = Ao/ except in an 
extremely small interval around r = 0. In that case, 

a = 26/(5+1) + 1/2 , (l»r»OA ) , (2.23c) 

which yields a w 2.17 if 5 w 5. 

Finally, at larger values of r, or sufficiently close to T c 
that -ff c is small enough to invalidate Eq. (|2.20(l . we are 
back to the paramagnetic case, Eq. I|2.19|l holds, and thus 
a = 3/2. 

One thus faces a rather complicated situation, where 
the exponent a can take on values between the Ginzburg- 
Landau value 3/2 and a value larger than 2, Eq. (|2.23c|l . 
depending on various parameters that are not easy to 
control or even determine experimentally. We will discuss 
this in more detail in Sec. IIVI Before we do so, in the 
following section we will give a more technical and more 
detailed derivation of all of our results. 



usual Ginzburg-Landau equations that includes the spin 
degrees of freedom. Far from magnetic criticality, the 
latter can be integrated out to yield ordinary Ginzburg- 
Landau theory with fi n entering the magnetic energy den- 
sity. At magnetic criticality, /i n becomes field dependent, 
which changes the thermodynamic critical field. In addi- 
tion, the leading term in the London equation vanishes, 
which leads to a generalized London equation that de- 
scribes exponential decay on a length scale shorter than 
Ao, in agreement with the qualitative arguments in Sec. 
Hll and with implications for the critical current as dis- 
cussed there. Unlike in the previous general discussion, 
in most of this section we will treat the magnetic critical 
behavior in a mean-field approximation. 



A. LGW theory for superconducting and magnetic 
fluctuations 



III. GENERALIZED GINZBURG-LANDAU 
THEORY 

We now consider a coupled field theory that describes 
both superconducting and spin degrees of freedom in 
order to derive the above results from a more micro- 
scopic level and gain a deeper understanding of their 
origin. Specifically, we consider a generalization of the 



Our starting point is an action for a complex scalar 
field ip describing the superconducting degrees of freedom 
coupled to a vector potential A, and a real vector field M 
describing the spin degrees of freedomii*ii We reiterate 
that the qualitative behavior near the superconducting T c 
does not depend on the symmetry of the order parameter, 
so our restriction to a scalar order parameter does not 
imply a loss of generality. The action reads 



S = 



dx 



1 

2m 
a 



(VM(i)) 



igA(*))V(x)| a + ||#iO| 



\^)\ A 



8tt 



B 2 (x) 



■M\x) 



"2 



y M 2 (x)) 2 - M(x) ■ B(x) - — H(x) ■ B(x) 



(3.1) 



Here and in the remainder of this section we use units such that Planck's constant and the speed of light are unity, 
H = c = 1. The first line is the standard Landau- Ginzburg- Wilson (LGW) functional for singlet superconductors. 
The first three terms in the second line are a standard vector-7W 4 theory, with M(x) the fluctuating magnetization. 
M couples to the vector potential via the M ■ B terrn^ with B = V x A, and the last term is necessary to relate S 
to the appropriate Gibbs free energy, see Eq. (|2.2() . Notice that ip and M are coupled only indirectly via the vector 
potential A. Spin-flip scattering of electrons by the magnetic moments does give rise to a direct coupling of the form 
M 2 \ip\ 2 £ but these terms are not important for our purposes. 

Minimizing this action with respect to ?/>*, A, and M yields the following saddle-point equations, 



h ip(x) + m ^(a;)! 2 ^^) - — (V - iqA(x)) 2 tp(x) 

m 

2 

-i — \ip*(x)Vil;(x) - i>(x)Vi>*(x)] - — \^(x)\ 2 A(x) 
2m m 



0. 



— V x \B(x) - H(x) - AttM(x)] , 

4-7T 



t 2 M(x)-aV 2 M(x)+u 2 M 2 (x)M{x) = V x A(x). 



(3.2a) 

(3.2b) 
(3.2c) 



r 



If we drop Eq. I|3.2c(l and put M = in Eq. (|3.2b|l superconducting solution of the full equations is tp = 0, 
(this corresponds to dropping M from the action) we B = H + 4ttM, and M determined by the magnetic 
recover the usual Ginzburg-Landau equations. 20 A non- 
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equation of state 

(r - aV 2 ) M{x) + u 2 M 2 {x)M{x) = H{x), (3.3) 

where r — ti — 4ir. For a small constant external field H 
a solution of Eq. I|3.3|l is M = \ n H, with 



(3.4) 



Xn = 1/r 



the normal-state magnetic susceptibility (see Ref. |25 
At this point it is the bare susceptibility, but it is clear 
that by renormalizing the spin part of the action before 
constructing the saddle-point solution one can make it 
the physical susceptibility. 



B. Effective theory for paramagnetic 
superconductors 

Now consider the full Eqs. (|3.2|) . For a small and slowly 
varying M(x) we have from Eq. (|3.2c|l 



M{x) = (x^ 1 + 47r) l B{x) 



(3.5) 



Substituting this into Eq. 13.2b|l we obtain 



j sc (a;) = V x B{x) - — V x H (x), (3.6a) 

47T/l n 47T 



where 



j sc (x) = -if- [f(x)V^(a!)-^x)Vf(x)] 

-— |^(a;)| 2 A(a3). (3.6b) 
m 

Together with Eq. (|3.2a(l , these are the equations of mo- 
tion for an effective action^ 



S, 



elf 



dx 



1 

2m 



(V -iqA(x))ip(x)\' 



t 



\^{xt 



\4>(x)\ 4 



_L_ B 2 (x) - — H(x) ■ B(x) 



(3.7) 



where we have dropped the now-superfluous subscript on 
the Landau parameters t and u. The same result is of 
course obtained by starting with Eq. (|3.1|l and integrating 
out M in a Gaussian approximation. 

The quantity j sc in Eqs. (|3.6|l is indeed the super- 
current, as can be seen by comparing Eq. Ij3.6a|) with 
Eq. (|2.9a|) . It does not explicitly depend on fi n , see Eq. 
(|3.6b[) . and this is important for the flux quantum to 
be independent of [i n . The magnetic energy _B 2 /87r/x n , 
which does explicitly depend on /i n , does not appreciably 
contribute to the free energy of a thin film or wire sample, 
and the standard determination of the critical current, 
Ref. |23, thus leads to the usual Ginzburg-Landau result 
with no correction due to /x n ^ 1. This corroborates the 
educated guess in Sec. Ill A 31 



For all other quantities, the usual analysis of Ginzburg- 
Landau theory now applies^ One characteristic length 
scale is given by the square root of the ratio of the coef- 
ficients of the gradient-squared term and the ip 2 term in 
Eq. H3.7f) . This is the superconducting coherence length 
£ = \Jl/m\t\. Another one is given by the square root of 
the ratio of the coefficients of the terms quadratic in A. 
For a constant ip, this is the London penetration depth 



A = y/m/iirq 2 ^ 2 ^ = Aq/VaV 



(3.8) 



This is identical with Eq. i|2.16bjl . which had been de- 
duced on elementary phcnomcnological grounds. 

For the Ginzburg-Landau parameter k = A/£ we now 
have k = kq/^/JI^, with kq the value of the parameter 
for ji n = 1. This implies that the superconductor is of 
type I or type II, respectively, for no < fx n /2 or ko > 
y/n n /2. While one can show this by an explicit analysis 
of the effective action, a fast way to relate the theory 
for arbitrary values of /i n to the one for /i n = 1 is to 
rewrite the action in terms of dimensionless quantities. 24 
In conventional Ginzburg-Landau theory, this is done by 
introducing 

x = X x , ip(x) = ipo 4>(x) , A(x) = V2H°\ A(x), . 

(3.9) 

Here ipo = y/—t/u is the superconducting order parame- 
ter scale. In terms of these quantities, the effective action 
reads 24 



— V - iA(x) I ip(x) 

K 



.\#@)\* + ±\fa)\* + ±. (v Ax, 

2 /i„ 



2H(x) • V x A(x) \ , 



(3.10) 



A simple further rescaling procedure shows that S e g de- 
pends only on a single dimensionless parameter, rather 
than the two parameters kq and /z n . Define 

x = x/n , A{x/kq) = A(x) , tp(x/K ) = ip(x). (3.11) 

Then 

2 



V - iA(x) \ ip(x) 



47T 

|V#)| 2 + ^hMi)| 4 + ^ (V Aix) 
(ji£H(x)) ■ ( 



-2-^= (^H(x)) ■ (V x A(x) 



(3.12) 



This shows that the theory with an arbitrary /i n maps 
onto ordinary Ginzburg-Landau theory with the replace- 
ments 



H 



/]j^.H. 



(3.13) 
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kq = \J [i n /2 thus indeed marks the demarcation between 
type I and type II superconductors, and the critical fields 
can be immediately obtained from the usual results at 
fj, n = For the thermodynamic critical field H c , the 
upper critical field H C 2, the lower critical field H c i, and 
the surface critical field H C 3 we obtain 



H C 3 



9(1*0/ JJQ H° 



ff(«o) 
1.695 H c2 . 



\/2k q 



(3.14a) 
(3.14b) 

(3.14c) 

(3.14d) 



where the universal function g has the limiting behavior 

jinx + 0.08 + 0(1 /x) for a;>l/v / 2, 
9( - X >-y for x = l/V2. 

(3.14b) 

If one neglects the weak dependence of g on its argument, 
H c i is approximately independent of /z n . 

C. Superconductors at magnetic criticality 

As one approaches a ferromagnetic instability, /i n keeps 
increasing and can no longer be treated as a constant. 
There are two effects that become important for our pur- 
poses. First, in a normal metal fi n becomes strongly field 
or induction dependent. At r = this dependence is non- 
analytic and described by the critical exponent S. Sec- 
ond, as r becomes on the order of £^/A (see Sec. IllBf) 
in a superconducting phase, the difference between /x n 
and fj, B can no longer be neglected. Related to this, the 
gradient squared term in Eq. i|3.2c|) must be taken into 
account. We now consider these effects, starting with the 
nonanalytic field dependence in the normal state. 

1. Thermodynamic critical field 

At magnetic criticality in the normal state, r — 0, one 
hasS 



X 



n (r = 0,H)= X o (h/H ) 



1/8-1 



(3.15) 



Here \o is a microscopic susceptibility, and Hq is a mi- 
croscopic field scale. M and, for small values of H , B are 



therefore proportional to H 1 / 5 , or H cx B s . For small B, 
the number fi n should thus be replaced by a function of 
B with the following leading i?-dependence, 



fin -» (H /B) 



5-1 



(3.16) 



with Hq = (47rxo) (5 / ( ' (5 ~ 1 - ) Ho- The magnetic energy cost 
of the flux expulsion that results from the formation of 
a Meissner phase (which equals minus the normal-state 
magnetic energy) is now obtained by using Eq. 
in Eq. ijTTjl . It is E m /V = HB/4tt - B s+1 /8ttH^ 

H^ 1/S H 1+1 ^/8tt. The condensation energy is still given 
by E con d/V = t 2 /4u, which yields 



n = r~ 

u 



S/(S+1) 



1 



H, 



| t |2W+i), ( 3.i7) 



The thermodynamical critical field is thus weaker than 
in the paramagnetic case, and the i-dependence is con- 
sistent with Eq. 12.21fl . By comparing with Eq. (|2.5|l . 
we see that with respect to the thermodynamical critical 
field, effectively scales like /z„ - l/|t| 2 C*-i)/(*+i) at 
magnetic criticality. 

Equations (|3.15l) through l|3.17|l hold also for small but 
nonzero values of r as long as one is in the field scaling 
regime, i.e, as long as H in appropriate units is large 
compared to r to an appropriate power. We will discuss 
this in more detail in Sec. IIVI At this point we only 
mention that, since H c vanishes as \t\ — > 0, sufficiently 
close to T c one will lose the field scaling for any nonzero 
value of r, and H c will be given by Eq. (|3.14a|l . 



2. Generalized London equation 

The ordinary London equation is obtained from Eq. 
(|3.2bl) by dropping M(x) and treating ip(x) = i\> as a 
constant (London approximation). With V x = 

this leads to 



-Aq 2 B(x) = V x V x B(x). 



(3.18) 



Now take M into account. Using Eq. I|3.2c|) in Eq. (|3.2bp . 
we can eliminate B and derive an equation for M . Once 
M is known, B follows from Eq. (|3.2c|l . Within the 
London approximation one finds 



M(x) = -(A^ n )VxVxM(i)+^) V 2 M(a;) + (|^) A 2 V x V x V 2 M(x) - u M 2 {x)M(x) 

-uAqV x V x M 2 (x)M(x). (3.19) 

I 

Here fi n = (47r+r)/r as in Sec. lIII Bl ^ = \/4ir + r = As long as /i n « 1, the first term on the right-hand side 
y / a/(4n + r), and u = u/(iir + r). of Eq. (|3.19|) leads to a variation of M on a length scale 
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A = Xo/^/JI^. The second term is a small correction to 
the first one since ^ < Ao. So is the third term, which is 

of order (j^A V 2 ~ /A 2 , < 1 relative to the first 

one. The linearized version of Eq. I|3.19[) thus reduces 
to the ordinary London equation, Eq. I|3.18|) . with Ao — > 
A. However, for r = the first term vanishes. This 
makes the second term the leading one, and the third 
term, which is of order AqV 2 compared to the second 
one, cannot be neglected either. The linearized equation 
thus reads 

M(x)= ((S,) 2 V 2 [l + ^VxVx]M(x). (3.20) 

With the same interface geometry as in Sec. HEP this 
takes the form 

M(x) = (a)V(z)- (&)\mW(x). (3.21) 

This linear quartic ODE is solved by an exponential 
ansatz, M(x) = M^e^ px . The real solution that falls 
off for x — > oo shows damped oscillatory behavior. From 
Eq. (|3.2c|l we see that B(x) shows the same behavior as 
M(x), up to corrections of 0(£^/Ao). With the bound- 
ary condition B(x = 0) = fi n H we finally obtain 



B{x) = p n He 



cos 



(3.22) 



This is the solution of the linearized version of Eq. 
(|3.19() at r = 0. In addition to leaving out the terms of 
0(M 3 ), we have also ignored the fact that the permeabil- 
ity, whether p n or /j, a , does depend on B or M at magnetic 
criticality. In a mean-field approximation, /i n cx 1/B 2 at 
r = 0, see Eq. i|3.16|l . which also leads to terms of 0(A1 3 ) 
in the nonlinear equation. Depending on the ratio of the 
external field to Hq, these terms may or may not be im- 
portant for the initial decay of M or B near the normal 
metal-to-superconductor boundary. However, once M or 
B has decayed sufficiently, these terms become sublead- 
ing compared to the linear ones in Eq. (|3.21|) . and the 
asymptotic behavior as B — > is always given by Eq. 

In order to make contact with the discussion in Sec. 
Ill Bl for small but nonzero values of r, let us consider the 
linearized Eq. 13.19fl while keeping the first term. Instead 
of Eq. (|3.21|l we then have 

M(x) = {xl/n n + (e? n ) 2 ) M"{x) - (&)\m™(x). 

(3.23) 

This is solved by 

M(x) = M e px , (3.24a) 

with 



P 2 = 



2A 2 , e° 



-4A 



(^m) 



(3.24b) 



Here we have chosen the solution for p 2 that yields 
p 2 — » 1/Ag for r — > oo. Equation l|3.24b|l still provides 
two solutions for p, and the physical solution for M is 
determined by the requirement that M be real. 

A discussion of Eq. Ij3.24b() shows that p 2 becomes 
purely real and negative at r = r s = + 

O ((ji^j /Ao^j • This is in agreement with the results of 

Blount and Varma, » who showed that spiral magnetic or- 
der coexisting with the superconductivity occurs at this 
point. For \r\ <C Cm/^o one has p 2 w — i/^Xo, which 
leads to Eq. IpT^ . For r > i°JX one finds p 2 » p n /Xl, 
which leads to 

B(x) = n n He- x ^ Xo , (3.25) 
in agreement with Eq. (|3.8|) . 

3. Penetration depth, and critical current 

Equation (|3.22|) shows that the effective penetration 
depth at magnetic criticality is 



A = V2£°A , (|r|«OAo), (3.26) 

in agreement with the conclusions of Ref. drawn from 
studying the ferromagnetic phase, and with Eq. (|2.22() 
with 7=1. The latter approximation results from the 
fact that our saddle-point equations of motion describe 
the magnetic equation of state in a mean-field approxima- 
tion. The discussion of Eq. (|3.24b|) shows that this result 
is valid for \r\ -C ^/Ao- By comparing with Eq. (|2.16b|) 
or H3.8|) . we see that with respect to the penetration 
depth, p, n at magnetic criticality scales like fi n ~ 1/ y'jtf 
in mean-field approximation, or p n ~ l/|i| 7 / 2 in general. 
The fact that 1/yT^n in Eqs. (|2.5(l and l|2.16b|l . respec- 
tively, must be interpreted differently for fj, a — > oo is a 
consequence of the influence of the superconductivity on 
the spin response. 

For r > £°,/Ao we have, from Eq. l(3~2T)| 



A — Ao/ 



(3.27) 



in agreement with Eq. (j3.8|l . 

The expression for the critical current given by Eq. 
(|2.18l) is general within the London approximation. Wc 
have now given a derivation of the behavior of the ther- 
modynamical critical field and the penetration depth 
given on phenomenological grounds in Eqs. (|2.21|) and 
(|2.22|) . respectively. The behavior of the critical current 
at or near magnetic criticality is thus given by Eqs. H2.23fl . 



4- Critical field H c z 

The critical exponent 7 is positive (7 ~ 1.4 for typical 
ferromagnetic universality classes in three dimensions 2 ^) . 
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The result for A, Eq. (|2~2l|) or lj3~2"3"j) in mean-field ap- 
proximation, of the previous subsection therefore means 
that A diverges more slowly for \t\ — » than the super- 
conducting coherence length £ oc l/y/\t\. Consequently, 
superconductors at magnetic criticality (|r| <§C Cm/^o) 
are necessarily of type l£ 

This observation notwithstanding, the critical field 
H c2 , which in a type-II superconductor signalizes the 
boundary of the vortex phase, still has a physical mean- 
ing: It is the minimum field to which the normal metal 
can be 'supercooled' before it discontinuously develops 
a nonzero superconducting order parameter S It is thus 
still of interest to determine i? C 2- Furthermore, the be- 
havior will be necessarily of type I only for \r\ in an ex- 
tremely narrow region. Outside of this region, Eq. (|3.27|) 
holds, and for a sufficiently large value of k = Ao/£ the 
superconductor will still be of type II. The determination 
of H C 2 is done by linearizing the Ginzburg-Landau equa- 
tion, Eq. (|3.2a[l . in ip. It then turns into a Schrodinger 
equation for a particle in a vector potential A, with 
—ti/2 = — 1/2 playing the role of the energy eigenvalue. 
By means of standard arguments^ this leads to a critical 
value of the magnetic induction B — V x H given by 
B C 2 = H® 2 = —tm/q. In a paramagnetic superconductor, 
this leads to 

Hc2 = H° 2 /fi a {fx n = const.), (3.28) 

which is the same as Eq. (|3.14bp . At magnetic criticality, 
we have, cf. Eq. (|3.1t)|) . 

H c2 = B & c2 /H^- 1] cx \t\ & . (3.29) 

Notice that, in this context, /i n scales as /i n ~ l/j^ 15-1 , 
whereas it scales as fi n ~ l/\t\ if the relevant field scale is 
H c . Since H C 2 vanishes much faster than H c , Eq. I|3.17[) . 
the field scaling region will be restricted to larger values 
of \t\, and H c2 will be given by Eq. (|3.28|l in a substantial 
range of t- values. Will come back to this in Sec. IIVI 

IV. DISCUSSION AND CONCLUSION 

To summarize, we have determined the electrodynamic 
properties of superconductors close to a ferromagnetic 
instability, i.e., materials that, in the absence of super- 
conductivity, would be paramagnetic with large ferro- 
magnetic fluctuations. This work complements previous 
studies of the coexistence of superconductivity with fer- 
romagnetic order»i*2 We have treated the superconduc- 
tivity in mean-field (Ginzburg-Landau) approximation. 
In addition, we have employed the London approxima- 
tion, treating the superconducting order parameter as 
a constant. The ferromagnetic critical point we have 
treated explicitly in a mean-field approximation, and we 
have used scaling arguments to consider the consequences 
of the exact magnetic critical behavior for the super- 
conductivity. We have found that the thermodynami- 
cal critical field H c decreases due to the ferromagnetic 



fluctuations, as one would expect, and depends on the 
magnetic critical exponent 5, see Eqs. (13.17(1 and l(2.21|l . 
However, the London penetration depth also decreases, 
which is intuitively less obvious. At magnetic criticality 
the behavior of the magnetic induction at a vacuum-to- 
superconductor (or normal metal-to-superconductor) in- 
terface is still characterized by exponential decay, but 
the characteristic length scale A is different from the 
usual London penetration depth Ao. Within a mean- field 
description of the magnetic criticality it is the geomet- 
ric mean of the zero-temperature magnetic correlation 
length and Ao, see Eqs. (|3.22|) and l|3.26[) : more gener- 
ally, it depends on the magnetic critical exponent 7, see 
Eq. (|2.22|1 . However, this behavior of the penetration 
depth is valid only within an extremely small region of 
width ^/Ao around magnetic criticality. Outside of this 
region, but still within the ferromagnetic critical region, 
the temperature dependence of the penetration depth is 
the same as in Ginzburg-Landau theory, see Eq. H3.27fl . 
For the critical current j c cx H c /X this implies a depen- 
dence on the reduced temperature given by \t\ a , where 
the exponent a depends on both 5 and 7, or on 5 only, 
depeding on the value of r, see Eqs. (|2.23|) . With expo- 
nent values appropriate for the usual ferromagnetic uni- 
versality classes, a w 1.8 extremely close to magnetic 
criticality, and a ~ 2.15 somewhat farther away. 

Let us now discuss these results in some more de- 
tail, and relate them to the experimental observations 
reported in Ref . lift 

For the temperature dependencies of various observ- 
ables at magnetic criticality we have assumed that the 
system stays tuned to magnetic criticality while the tem- 
perature is varied. Let us discuss to what extent this 
assumption is realistic. Consider a phase diagram in a 
plane spanned by the temperature and some non-thermal 
control parameter x, e.g., the hole doping concentration 
in the case of MgCNi3 and consider the following two 
qualitatively different possibilities. Figure [21 shows a sit- 
uation where the magnetic phase separation line does not 
cross the line x = 0. The stoichiometric compound thus 
does not enter a magnetic phase upon cooling, although 
the system is close to a magnetic transition for all tem- 
peratures below the superconducting T c . This scenario is 
believed to apply to MgCNi 3 . Figure shows a situation 
where the magnetic phase separation line does cross the 
line x = 0, so that the stoichiometric compound enters 
a phase where superconductivity and magnetism coexist 
at some temperature below T c . This is the situation that 
was discussed in Refs. Q andy and observed in ErRli4B4 
and HoMoeSgi^ The magnetic transition is to a phase 
with spiral magnetic order at a temperature T s slightly 
below the temperature where ferromagnetism would 
occur in the absence of superconductivity. 1 

We now can see what is required to keep r constant 
while varying t, namely, a situation as shown in Fig. [21 
with the dashed line essentially parallel to the T-axis. r is 
then given by the dimensionless distance between the two 
lines. In order for the penetration depth to display the 
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T 

x = > 

T c 

FIG. 2: Schematic phase diagram showing a normal metal 
(NM), a ferromagnet (FM), a superconductor (SC), and a 
ferromagnetic superconductor (FMSC) in a temperature (T) 
- control parameter (x plane. The solid line denotes the su- 
perconducting transition, the dashed line, the magnetic one. 
Along x — there is only one phase transition at the super- 
conducting T c . See the text for additional explanation. 




x = 1 ' i 

T s T m T c 

FIG. 3: Same as Fig. |21 but with a magnetic transition for 
x — at a temperature T m < T c . On the x = axis it is 
shown that T m splits into the bare magnetic transition tem- 
perature and the physical transition temperature T B to 
a state with spiral magnetic order, RefQ. See the text for 
additional explanation. 



non-Ginzburg-Landau behavior described by Eq. (|3.26l) 
or, more generally, Eq. (|2.22() . the two lines would have to 
be extremely close, in order to keep r smaller than ^/Ao, 
see Eq. <|3.26ll ■ This would result in a temperature depen- 
dence of the critical current given by Eqs. (|2.23all2.23b(l . 
While this is possible, it is a very non-generic situation, 
and it would result in a very large magnetic susceptibil- 
ity of the normal metal just above the superconducting 
transition temperature. 

A situation that is still very non-generic, but requires 
somewhat less fine-tuning, is one where the dashed line is 



still essentially parallel to the T-axis, but in a somewhat 
larger r-range, say, with r on the order of a few percent. 
In this case the penetration depth will show the usual 
1/lij 1 / 2 temperature dependence, see Eq. (|3.27|1 . The 
temperature dependence of the thermodynamic critical 
field will be more complicated in this case. The general- 
ization of Eq. (|2.20|l to nonzero values of r is 

Xn = r-^f x (H/r^), (4.1) 

with 7 = [3(8 — 1), [3, and 8 the usual critical expo- 
nents for the magnetic transition. In order for Eq. (|2.20|) 
to hold, the H must be large compared to r 136 in suit- 
able units. The latter are not determined by any uni- 
versal arguments, but an analysis of the critical equa- 
tion of state for both the high-temperature ferromag- 
net Ni (T m 630K)SS, and the low-temperature ferro- 
magnet CrBr 3 (T m « 33 K)2i shows that in either case 
the relevant energy or field scale (we use units such that 
A;b = fJ-B = 1) is given by T m , which is plausible. The 
crossover between the field scaling that leads to Eq. (|2.21() 
and the static scaling that leads to Eq. (|2.5|l thus occurs 
at a crossover field 

ffx«T°r^. (4.2) 

[35 ~ 5/3 for ferromagnetic phase transitions, and with 
T m w 10 K and r w 0.1, one finds H x w 0.02 T^. For 
MgCNi 3 in the vicinity of T c , this leads to H x w 0.2 T. 
With H C 2 at zero temperature on the order of 14 T and 
k w 40, 28 one expects H C (T = 0) = H c2 /V2k ps 0.25 T. 
Since H c vanishes at T c , this means that H c will be given 
by Eq. (|2.21|) sufficiently far away from T c , but cross over 
to H c oc \t\ near T c . Consequently, the critical current 
exponent a will be given by Eq. (|2.23c|) at some distance 
from T c , and cross over to the Ginzburg-Landau result 
a = 3/2 as \t\ — > 0. In the experiment of Ref. fl6L no such 
crossover was observed down to |i| « 0.01. 

At least within the London approximation, our results 
confirm the conclusion of Ref. |2| that superconductors 
near a ferromagnetic instability are necessarily of type 
I. However, we have also shown that this conclusion is 
inevitable only within an extremely small region around 
the (bare) magnetic critical point. The fact that MgCNi3 
is observed to be of type II2& is therefore not necessar- 
ily in contradiction to the notion that this material is 
almost ferromagnetic. However, Eq. I|3.29|l predicts a 
strong deviation from Ginzburg-Landau behavior for the 
upper critical field H c2 . Since H c2 goes to zero rapidly as 
\t\ — > 0, this behavior will show only at substantial values 
of \t\ even if r is very small. No anomalous behavior was 
observed for |i| up to 0.51^ This is reconcilable with close 
proximity to a magnetic instability only if r is very small 
close to T c , and grows with decreasing temperature, in 
which case H C 2 might never show the magnetic critical 
behavior. A signature of this situation would be a large 
magnetic susceptibility in the normal state just above T c . 

The conclusion from this discussion with respect to the 
experimental observations in Ref. ^| is as follows. While 
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it is possible that proximity to a ferromagnetic instability 
is the cause of the observed anomalous behavior of the 
critical current, such an explanation requires fine tuning 
of the phase diagram, and would have to be accompanied 
by a very large enhancement of the spin susceptibility in 
the normal phase just above T c . Explaining the lack of 
an anomaly in the temperature dependence of H C 2 prob- 
ably requires that the material is closer to the magnetic 
instability near T c than at T = (i.e., the dashed line 
in Fig. Incomes closer to the T-axis with increasing T). 
A direct measurement of the spin susceptibility in the 
normal phase would be of great interest in this context. 

Finally, we discuss our predictions for the case of a 
superconductor that does undergo a transition to a mag- 
netic state below T c , i.e., the situation represented by 
Fig. |3 In the (very small) temperature interval of width 
2|T^ — T s \ around T^, both the thermodynamic criti- 
cal field H c and the penetration depth A will show an 
anomalous temperature dependence, and the critical cur- 
rent exponent will be given by Eq. (|2.23b|) . Outside of 
this region, but not too close to T c , H c will be anoma- 



lous, but A will be conventional, and the critical current 
exponent will be given by Eq. (|2.23c[) . Upon approach- 
ing T c , H c will fall below the crossover field given by Eq. 
I)4.2|) . and its temperature dependence will cross over to 
the usual linear Ginzburg-Landau behavior. The critical 
current exponent close to T c will thus be the conven- 
tional a — 3/2. The location of this crossover depends 
on the critical field scale, and will thus be material de- 
pendent. Critical current measurements in the materials 
like ErRli4B4, or HoMogSs, which are believed to fall into 
this class, would be very interesting. 
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